SHY COUPLINGS 



Itai Benjamini, Krzysztof Burdzy and Zhen-Qing Chen 

Abstract. A pair (A, Y) of Markov processes is called a Markov coupling if A 
and Y have the same transition probabilities and (A, Y) is a Markov process. 
We say that a coupling is "shy" if there exists a (random) e > such that 
dist(At,l^) > e for all t > 0. We investigate whether shy couplings exist for 
several classes of Markov processes. 



1. Introduction. The proofs of the main theorems in two recent papers, [BC2] and 
[BCJ], contained arguments showing that certain processes come arbitrarily close to each 
other, at least from time to time, as time goes to infinity, with probability one. The proofs 
were based on ideas specific to the models and were rather tedious. We decided to examine 
several classes of Markov processes in order to determine the conditions under which there 
exists a pair of Markov processes defined on the same probability space such that each 
marginal process has the same transition probabilities and the two processes do not come 
close to each other at any time. Although we do not have a complete solution to this 
problem, we offer a number of results whose diversity points to a rich theory. Some of our 
theorems, examples and techniques may have interest of their own. 

We will focus on two classes of processes — refiected Brownian motions on Euclidean 
domains and Brownian motions on graphs. The second class of processes is really discrete 
in nature, in the sense that similar techniques work for random walks on graphs. We chose 
these classes of processes because similar processes appeared in our research in the past. 

For a general overview of coupling techniques, see [L] . 

The rest of the paper is organized as follows. We present basic definitions and elemen- 
tary examples in Section 2. Section 3 is devoted to Brownian motions on graphs. We show 
that there exists a shy coupling for Brownian motions on a graph if all its vertices have 
degree 3 or higher. Four examples are also given to illustrate the case when the graph has 
some vertices of degree one. Section 4 deals with refiected Brownian motions on Euclidean 
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domains, showing that there exist no shy couphngs on C"^ -smooth bounded strictly convex 
domains. 



2. Preliminaries and elementary examples. Unless specified otherwise, all pairs of 
processes (X, Y) considered in this paper will be "Markov couplings," i.e., they will satisfy 
the following assumptions. 

(i) {Xf, t > 0}, {y^, t > 0} and {{Xf, Yt),t > 0} are Markov, and the transition probabil- 
ities for X and Y are identical. 

(ii) The distribution of {Xt,t > s} conditional on {{Xs,Yg) = {x,y)} is the same as the 
distribution of {Xt, t > s} conditional on {Xg = x}, for all x, y and s. 

Our definition of a Markov coupling is slightly different from similar concepts in the 
literature. One could investigate the question of whether our results hold for "couplings" 
defined in other ways, for example, whether condition (ii) is essential. However, we feel 
that there are more exciting open problems in this area — see the end of Section 4. 

The following elementary discrete-time example shows that there exist couplings that 
satisfy (i) but do not satisfy (ii). 

Example 2.1. We take {0, 1} as the state space of a discrete time Markov process and 
we let {Xfe, A; > 0} be a sequence of i.i.d. random variables with P(Xfe = 0) = P(Xfc = 
1) — 1/2. We will define a process {Yk, A; > 0} with the same distribution as {X^, k > 0}. 
We let Yq be independent of {X^, k > 0}. For /c > 1, we construct Yk so that P(yfc = 
I Xk-i = 0) = 0.7 and P(Ffc = 1 | Xk-i = 1) = 0.7. Moreover, for every A; > 1, we 
make Y^ independent of X/s ior j < k — 1. It is elementary to check that {X^, > 0}, 
{^A:, ^ > 0} and {(Xfc, Y^), k > 0} are Markov but for j > 1, the distribution of {Yfc, k > j} 
conditional on [Yj-i = 0} is not the same as the distribution of {Yk,k > j} conditional 



We will assume that the state space S for Markov processes X and Y is metric and 
we will let d denote the metric. The open ball with center x and radius r will be denoted 
B{x,r). The shortest path between two points in S will be called a geodesic. For some 
pairs of points, there may be more than one geodesic joining them. 

Definition 2.2. A coupling (X, Y) will be called shy if one can find two distinct points x 
and y in the state space with 



on{(X,_i,Y,_i) = (0,0)}. 



P inf d{Xt, Ft) > I Xo = X, Fo = y > 0. 
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Note that the term "shy couphng" is a label for a family of Markov transition proba- 
bilities. 

We proceed with completely elementary examples of shy and non-shy couplings. 

Examples 2.3. (i) Let X be a Brownian motion in and let 7^ y e R'' be a fixed 
vector. Let Yt = Xt + y for all t >0. Then (X, Y) is a shy coupling. 

(ii) Let X be a Brownian motion on the unit circle in and let 6 e (0, 27r) be a fixed 
number. We define Yt using complex notation, Y^ = e^^X^ for all t > 0. Then (X, Y) is a 
shy coupling. 

(iii) The last two examples can be easily generalized to a wide class of Markov processes 
on spaces S with a group structure. If there is a group element a ^ such that a + X is 
a Markov process having the same transition probabilities as X and inffe^^ d(6, a + b) > 
then (X, a + X) is a shy coupling. 

(iv) Let X and Y be independent Brownian motions in R^. Then (X, y) is a shy 
coupling if and only if d > 3. 

In the sequel, for a, 6 G M, a A 6 := min{a, 6}, a V 6 := max{a, b} and a"*" := a V 0. For 
a > 0, [a] denotes the largest integer that does not exceed a. 

3. Brownian motion on graphs. In this section, wc will consider processes whose state 
space is a finite or infinite graph. More precisely, let Q = (V, S) be a graph, where V is the 
set of vertices and S is the set of edges. We will assume that all vertices have a finite degree, 
i.e., for every vertex there are only a finite number of edges emanating from this vertex, 
but we do not assume that this number is bounded over the set of all vertices. Wc allow 
an edge to have both endpoints attached to one vertex. Every vertex will be attached to 
at least one edge. We will identify edges with finite open line segments (connected subsets 
of R) , with finite and strictly positive length, and we will identify vertices with topological 
endpoints of edges. In this way, we can identify the graph Q with a metric space [S, d), 
where S = S UV, and d{x, y) is the shortest path between x and y along the edges of the 
graph. We will assume that the length of any edge is bounded below by ro > 0. 

Next we will construct "Brownian motion" X on S. See [FW] for a definition of a 
general diffusion on a graph. We leave it to the reader to check that our somewhat informal 
description of the process is consistent with the rigorous construction given in [FW]. By 
assumption, our process will be strong Markov. Suppose that x E e E £ and x is not an 
endpoint of e. Recall that e can be identified with a line segment, say, e = [0,y]. Then 
X e (0, J/). If Xq = X, then the process X evolves just like the standard one-dimensional 
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Brownian motion until the exit time from (0,j/). Next suppose that a; e <S is a vertex. 
Then there are n edges ei, 62, . . . , e^, attached to x, with n > 1. Choose a small r > such 
that the ball B{x, r) consists of line segments Ij, j = 1, 2, . . . , A;, which are disjoint except 
that they have one common endpoint x. Note that n < k < 2n, but not necessarily k = n, 
because some edges may have both endpoints at x. We will describe the evolution of X 
starting from x until its exit time from B{x, r). Generate a reflected Brownian motion R 
on [0, 00), starting from 0, and kill it at the first exit time from [0, r], denoted T^. Label its 
excursions from with numbers 1, 2, . . . , /c, in such a way that every excursion has a label 
chosen uniformly from {1, • • • , fc} and independently of all other labels. Then we define 
Xt for t e [0, Tr] so that d{Xt, x) = Rt and Xt G /j, where j is the label of the excursion 
of R from that straddles t (if Rt = then obviously Xt = x). This defines the process 
Xt until its exit time from B{x,r). What we said so far and the strong Markov property 
uniquely define the distribution of X . Note that when the degree of a vertex x is 1 then 
X is best described as a process refiected at x. The process X spends zero amount of time 
at any vertex. 

Recall that we have assumed that the length of all edges is bounded below by tq > 0. 
Under this assumption, the process cannot visit an infinite number of vertices in a finite 
amount of time. Hence, the above construction defines a process for all t > 0. Another 
consequence of the assumptions that all edges have length greater then tq and all vertices 
have finite degree is that for any two points in S there is only a finite number of geodesies 
joining them. It is clear from our construction that vertices of degree 2 will play no essential 
role in the paper and can be ignored. So we will assume without loss of generality that 
there are no vertices of degree 2. 

Theorem 3.1. If all vertices of Q have degree 3 or higher then there exists a shy coupling 
for Brownian motions on S. 

Proof. We will construct a coupling (X, Y) of Brownian motions on S such that X and Y 
move in an independent way when they are far apart and they move in a "synchronous" way 
when they are close together. Clearly, independent processes do not form a shy coupling on 
a finite graph. Remark 3.2 below explains why it is hard, perhaps impossible, to construct 
a "synchronous" shy coupling. 

For any x,y E S with d{x, y) > ro/4, we will define {Xt, Yt) starting from (Xq, Yq) = 
{x,y), for t e [0,t], where r is a random time depending on x and y. Then we will 
explain how one can define {Xt,Yt) for t e [0, 00) by pasting together different pieces of 
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the trajectory. 

(i) Recall that the length of any edge is at least ro > 0. First suppose that d{x, y) > 
3ro/4. Then we let {{Xi^Yi)^t e [0,t]} be two independent copies of Brownian motion on 
S and we let r = inf{t > : d{Xt,Yt) = ro/2}. 

(ii) Next suppose that x,y e S are such that d{x,y) G (ro/4, 3ro/4), and none of these 
points is a vertex. Let 

a(r)= ^^'^'-^°)^ Al, (3.1) 

and B and B' be independent Brownian motions on R starting from the origin. Let Ut = Bt 
and 

dVt = ^/l-a''{Ut-Vt)dBt + a{Ut - Vt)dB'„ (3.2) 
with Vq = vo = d{x,y) > ro/4:. Then, if we write Zt = Vt — Ut, we obtain 

Zt = Vo + [ i^/l-(7^{Zs)-l)dBs+ I (T{Zs)dB'„ 
Jo Jo 

and for Z^ = Zt - ro /4, 

Z^^vo-ro/4+ [ (Vl-cT2(Z^ + ro/4)-l)dS,+ [ a{Z'^ + r^/ 4)dB'^. 
Jo Jo 

So 



Zi^vo-ro/4+ [ 'y{Z'^)dWs 
Jo 



where 

df 



7(r) = (Vl-^'(r + ro/4) - 1)^ + a\r + ro/4) 

and is a Brownian motion on M with l^o = 0. The process Z' has the same distribution 
as 

vo -ro/4 + Wrt, 

where 

Tt:=mi\s>0: / 'yivo - ro/4 + Ws)~^ds > t 



Note that for small r > 0, 7(r) = O(r^) and so in particular J^^ 7(r) ^dr = oo. Thus by 
Lemma V.5.2 of [KS], 

/ 7(^0 - ro/4 + Ws)~^ds = oo 

^0 



almost surely, where Tq = inf{t > : vq — ro/4 + = 0}. We conclude that Z' never hits 
0; in other words, Z never reaches ro/4. 

Suppose that Xq = x, Yq = y and recall that we have assumed that d{x,y) e 
(ro/4, 3ro/4). Suppose that x E ei E £ and y E €2 E £. Let ei \ {x} consist of two 
line segments ef and el, with being the one closer to y. Similarly, 62 \ {y} consists of 
two line segments 62 and 62, and 62 is closer to x. We will define X and y on an interval 
[0, r] to be specified later. We define Xt on ei to be such that d{Xt, x) = \Ut\ and Xt e 
if and only if Ut < 0. We define the process Yt on 62 by conditions d(Yt, y) = | — vo \ and 
Ft e 63 if and only if Vt < vq. We let r be the first time t > that Xt or is at a vertex, 
or d{Xt, Yt) = 3ro/4. We see that over the interval [0,t), the distance between Xt and Yt 
remains in the interval (ro/4, 3ro/4). 

(iii) This part of our argument is based on the "skew Brownian motion." The skew 
Brownian motion U is a real-valued diffusion which satisfies the stochastic differential 
equation 



where S is a given Brownian motion with Bq = 0, (3 E [—1, 1] is a fixed constant and L 
is the symmetric local time of U at 0, i.e.. 



The existence and uniqueness of a strong solution to (3.3)-(3.4) was proved in [HS]. In the 
special case of = 1, the solution to (3.3) is the reflected Brownian motion. An alternative 
way to define the skew Brownian motion is the following. Consider the case (3 > 0. Take 
a standard Brownian motion B[ and flip every excursion of B[ from to the positive side 
with probability P, independent of what happens to other excursions (if an excursion is on 
the positive side, it remains unchanged). The resulting process has the same distribution 
as U defined by (3.3)-(3.4). For more information and references, see recent papers on 
skew Brownian motion, [BCl] and [BK]. 

Suppose that x,y E S, d{x,y) G (ro/4, 3ro/4), and a; is a vertex. Note that y is 
not a vertex. By assumption, the degree k of vertex a; is 3 or greater. Suppose that 
Bt is a Brownian motion on R and let U he a solution to (3.3)- (3.4), with /? defined by 
(1 — /?)/(! + /?) = A; — 1. Note that /? < 0. We label negative excursions of U from with 
numbers 1,2,..., A; — 1, in such a way that every excursion has a label chosen uniformly 
from this set and independently of all other labels. 



Ut = Bt + /3L: 



(3.3) 




(3.4) 
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Suppose that B' is a Brownian motion independent of B. Recall the definition of the 
function a and the process V given in (3.1) and (3.2), respectively, with Vq = vq = d{x,y). 
Then, if we write Zt = Vt — Ut, we obtain 

Zt = vo-PL^+ [ {^yi-a^iZs)-l)dBs+ faiZs)dB'^, 
Jo Jo 

and for = - ro/4, 

Zi = vo-ro/^-PL^ + [ {^/l-(7^{Zi + ro/^)-l)dBs+ f (j{Z'^ + ro/A)dB'^. 

Jo Jo 

We have already pointed out that for small r > 0, 

7(r) = (Vl-^'(r + ro/4) - 1)^ + ^^(r + ro/4) = ©(r^). 

Since /? < 0, the process —fiL^ is nondecreasing. These observations and the argument 
used in the first half of (ii) imply that Z' never hits 0, i.e., Z never reaches ro/4. 

The ball 3ro/4) consists of line segments /j, j = 1, 2, . . . , k. We assume that Ik is 
the line segment containing y. We define Xt on these line segments so that d(Xt, x) = \Ut\. 
If Ut > then Xt E Ik- If Ut < then Xt G Ij, where j is the label of the excursion of 
U straddling t. Suppose that y E e E £. Let e \ {y} consist of two line segments and 
e^, with being the one closer to x. We define Yj on e by d(lt, y) = \Vt — vo\. We let 
Yt E if and only if Vt < vq. We let r be the infimum of t such that Yt is at a vertex, or 
d{Xt,Yt) — 3ro/4. Observe that over the interval [0,t), the distance between Xt and Yt 
remains in the interval (ro/4, 3ro/4). 

Now we will define the process {Xt,Yt) for all t > 0, assuming that Xq = x, Yq = y 
and d{x,y) > ro/4. We use one of the parts (i)-(iii) of the proof to define the process 
(X, Y) on an interval [0, ti]. Then we proceed by induction. Suppose that the process has 
been defined on an interval [0, Tfe] and d(XT-^,yT-^) > ro/4. Then we use the appropriate 
part (i)-(iii) of the proof to extend the process, using the strong Markov property at r^, to 
an interval [0,rfc+i]. It is easy to see that — > oo a.s., so the process {Xt, Yt) is defined 
for all t > 0. It is straightforward to check that {Xt, t > 0} and {Yt, t > 0} are Brownian 
motions on S and {X, Y) is a shy Markov coupling, as defined in Section 2. □ 

Remark 3.2. One may wonder whether it is possible to construct the shy coupling in the 
proof of Theorem 3.1 using the skew Brownian motion in such a way that the distance 
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between X and Y does not change on time intervals where both processes stay away from 
vertices ("synchronous couphng"). This idea works well for many graphs but runs into 
technical problems when we have a configuration similar to that in Fig. 3.1, with many 
geodesies joining two vertices. Suppose that the lengths of edges in Fig. 3.1 are chosen so 
that there are 6 geodesies between x and y. We will argue that if Xq = x and Yq = y 
then for small t > the distance between X and Y has to decrease. This is because X 
will have to move towards z with probability 5/7, and Y will have to move towards z 
with probability 1/3. Since 5/7 + 1/3 > 1, X and Y will find themselves on a geodesic 
from X to y, moving away from their starting points towards each other, with positive 
probability. For this reason, we could not find a "synchronous" shy coupling based on the 
skew Brownian motion. 




Figure 3.1. 



The rest of this section is devoted to graphs that have at least one vertex of degree 
1 (recall that vertices of degree 2 can be ignored and we assume that Q does not contain 
any of them). We do not have a general theorem covering all graphs with some vertices of 
degree 1 but we have four examples illustrating some special cases. 

Example 3.3. This example is similar to Examples 2.3 (i)-(iii). Suppose that there exists 
an isometry I : S ^ S such that infj;^^ d(x, /(x)) > 0. It is not hard to show that this 
holds if the isometry has no fixed points, i.e., if there does not exist x & S with I{x) = x. If 
such an isometry exists, then we can first construct the process X and then take Yt = I{Xt) 
for all t > 0. Obviously, thus constructed coupling {X,Y) is shy. Fig. 3.2 shows that a 
graph with some vertices of degree 1 may have this property. 
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Figure 3.2. 



Our next lemma is a large deviations-type estimate. Recall that all edges are at least 
ro > units long, by assumption. 

Let X be Brownian motion on a graph S. For A G S, define Ta = inf {t > : Xt E A}. 

Lemma 3.4. Assume that the degrees of all vertices are bounded above by mo. There 
exist constants Cq > 0, to < oo depending on tq only such that for t G (0, to) and r > 
with r'^ > t and B{x, r)'^ 7^ 0, we have 



In applications of the above estimates, r will be "fixed" and then a t (much smaller 
than r^) will be chosen. For this reason, we did not try to optimize the non-exponential 
factors — most likely they are not best possible. 

Proof, (i) First we will prove the lower bound. We start with some preliminary estimates. 

Suppose that S is a Brownian motion on M with Bq = and let = inf{t > : 
Bt = r}. Then P(T^ < t) = 2^^= e~" /"^^du for r > 0. The following inequalities are 
well known (see Problem 9.22 on page 112 of [KS]): 

" -'^V2< re--'/'du<^-e-^'/'. 



l + r2 

So for r > 1, 



2r J J. r 
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By scaling we obtain for t <r'^, 



t 1 




2iT r 



exp 



r 

2i 



< P{T^ < t) 



oo 



r/Vt 

■2 /n , df 



e-'' l^dv < \l exp - — 

TT r \ 2t 



(3.5) 



For r < 1, we have a trivial lower bound e~" ^'^du = cq > 0. For t > r"^, the same upper 
bound holds but the lower bound has to be replaced with a trivial bound ci = \J~^cq > 0. 
Let to < oo be the largest real such that 

l-2exp(-r2/(2to)) > 1/2. 



We will derive an estimate for P(T^ < t AT 



-ro/2) 



for r > ro/2 and t < to. If < to, 



then r < C2ro for some constant C2 < oo. In this case, it follows easily from the support 
theorem for Brownian motion that P(T^ < t A T^^^^^2) > C3 > for every t G (r^, to). 

Now suppose that t < A to. If Brownian motion hits —ro/2 and then it reaches r 
in t seconds or less, it has to go from level —ro/2 to level r in t seconds or less. Hence, by 
the strong Markov property applied at T^^^j^, 



P{T^ < t ATf, 



•ro/2. 




>p(T,^<t)-p(r,^+,„/2<t) 

[2t 1 ( (r + ro/2)2\ 




(3.6) 



Suppose that yi is a vertex, d(?/i,?/2) = '''0/2, and lies between and j/2 so that 
d(yi, 2/3) + d(y3, ^2) = ro/2. Let d(yi, ^/3) = ri e [0, ro/2] and suppose that Xo = yz- The 
process i?t = d(Xt,yi) is a one-dimensional reflected Brownian motion with Ro = ri, at 
least until it reaches ro. Let {zi, Z2, ■ ■ ■ , z^} be the set of all points with d{zj,yi) = ro/2 
(j/2 is one of these points). Let X^/2 = inf{t > : -Rt = ro/2}. If X visited j/i before 
then it is at y2 at time X^/2 with probability 1/k, by symmetry. Since k < mo, the 
probability that X starts at y^ and reaches y2 in s seconds or less is greater than or equal 
to the probability that reflected Brownian motion that starts from ri reaches ro/2 in s 
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seconds or less, divided by mo- The last probability is bounded below by the analogous 
probability for the non-reflected Brownian motion, so using (3.5) we obtain, 

P(r- < . I X„ = ^3) > -,/f exp (Jl^lllpl) , (3.7) 

mo V 27r ro/2 — ri \ 2s J 

for s < {ro/2 — ri)^. If s > (ro/2 — ri)^, the bound is c\/mQ. Note that these estimates 
hold for all ri e [0,ro/2], including ri = 0. 

Consider any x' e dB{x, r) and let F C 5 be a geodesic connecting x and x' . Suppose 
that r contains some vertices and denote them a;i,X2, . . . -.Xk-, in order in which they lie 
on r, going from x to x' . Let xq = x and Xk+i = x' . If there is a vertex closer to Xq than 
ro/2 and it is not Xi then we let uq be the point at the distance ro/2 from that vertex, 
between xq and xi. For every xj, j > 1, we let yj G F be the point ro/2 away from xj, 
between Xj and Xj+i. 

Let J = 1, . . . , mi, be the sequence of all points Xj and yj, in the order in which 
they appear on F from x to x' , including x and x' . Note that ^i<j<^^_i d(2:j, zj+i) = r. 
By the strong Markov property applied at the hitting times of ^j's, the probability that 
X starting from x will hit x' in t seconds or less is bounded below by Hjl^ ^ Pj^ where 
Pj is the probability that X starting from Zj will hit Zjj^i in tj seconds or less, and 
tj = td{zj,Zj+i)/r. 

If Zj is a vertex or xq then, by (3.7), if tj < d{zj, Zj+i)'^, 

- mo V 27rd(;2„^,+i) ^ \ 2tj J ^ ^ 

If tj > d{zj, Zj+i)'^ then 

Pj > ci/mo > (ci/mo) exp {-^^^^^^) ■ (3-9) 
For other Zj^s we use (3.6) to obtain, for tj < d{zj, Zj+i)'^ , 



and for tj > d{zj, Zj+i) 



2 



p,>C3>C3exp(-^^|±ll!'). (3.11) 
11 



The product of exponential factors on the right hand sides of (3.8)-(3.11) is equal to 
n'-p('-%?±i£)=exp(-^). (3.12) 



2ti J \ 2t, 

If > d(2;j, 2:^-1-1)^ then the non-exponential factor in (3.9) is ci/rriQ and it is C3 in 
(3.11). The non-exponential factors in (3.8) and (3.10) are bounded below by 



1 jt, 1 _ 1 Jtd{z,,zj+i)/r 1 Ul_ 



mo V 27r 2d(2j, 2j_|_i) mo V 27r 2d{zj^ Zj-^i) mo V 27r 2r 

We conclude that the product of non-exponential factors in (3.8)-(3.11) is bounded below 

by 

I \ [r/ro] 

C4 t 1 




mo V 2tt 2r 

This combined with (3.12) gives the lower bound in the lemma. 

(ii) Next we will prove the upper bound. Let {r^} be the family of all Jordan arcs in 
S linking x with dB{x,r). The number of edges in B{x,r) is bounded by m2 = mo^^°^ so 
the number of Fj's is bounded by m^ = m2!. The length of any Tj is at least r. 

Consider some F^. We will define a process that measures the distance from Xt 
to X along Ffc, in a sense. We will "erase" excursions away from Ffc and loops as follows. 
For t > 0, let i{t) = sup{s < t : Xg E Tk}. Let Vfc = V \ F;;;, i.e., Vfe is the set of vertices 
that do not belong to F^. Let Ti = 0, 

Sj = mi{t > Tj : Xt e Vfc), j > 1, 
T,+i=inf{t>5,:X, = X,(s.)}, j>l. 

If t e [S'j,Tj_|_i] for some j > 1, we let be the distance from X£(^Sj) to x along Ffc. For 
other t, we let R^ be the distance from -^.g(t) to x along F^. The process R^ is a time- 
change of reflected Brownian motion, that is, it is reflected Brownian motion "frozen" on 
time intervals when X is outside F^ (and some other intervals). Hence, the probability 
that Rf reaches r in t seconds or less is less than the right hand side of (3.5). Note that 
one of the processes R^ must be at the level r at the time when X hits dB{x,r). Hence, 
an upper bound on the probability in the statement of the lemma is the product of the 
right hand side of (3.5) and m^. □ 

Example 3.5. Suppose that the graph S is compact and has the following structure. For 
some X E S, the set S \ {x} is disconnected and consists of a finite number of disjoint finite 
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trees 7i, 72, . . . , 7^, and a graph lA (not necessarily a tree). We say that a vertex of a tree 
is a leaf if it has degree 1. Assume that for some ri > r2 > and every leaf y ^ x oi any 
tree Tj we have d{x, y) > ri, and for any z eU, d{x, z) < r2 (see, for example. Fig. 3.3). 
Suppose that (X, Y) is a coupling of Brownian motions on S. We will show that (X, Y) is 
not a shy coupling. 




Figure 3.3. 



(i) Let T = [jj Tj \ {x}. First, we will show that there exist pi > and a stopping 
time Ti < oo such that with probability greater than pi, either Xt-^ = x and ^ ^7 
or = X and Xt^ € T. Let rs G (r2,ri) and W = {y E T : d(y, x) > r^}. It is easy 
to see that X is recurrent so T2 = inf{t > : Xf E W} < 00 a.s. Suppose first that 
Yt2 E T U {x} and let T3 = inf{t > T2 : Xt = x or Yt = x}. Then Ti = T3 has the 
properties stated above. 

Next suppose that Yt^ E U and let T4 = inf{t > T2 : Xt = x} and T5 = inf{t > T2 : 
Yt = x}. By Lemma 3.4, for some p2, > and all y G W and z E W, 

P(T5 < T2 + to I Xt^ = z, Yt, =y)> P(T4 < T2 + to I Xt, = z, Yt, = y) +P2. 

Hence, 'P{T^ < T4 \ Xt, = z,Yt, = y) > P2^ and it follows that we can take Ti = T5 on 
the event {T5 < T4}. This completes the proof of our claim, with pi = P2- 

(ii) Recall that the length of any edge is bounded below by tq > 0. Fix an arbitrarily 
small e E (0, ro/3). We will show in the remaining part of the proof that X and Y come 
within e distance to each other in finite time almost surely, which will then imply that 
(X, Y) is not a shy coupling. 

The rest of the proof is based on an inductive argument. We will now formulate and 
prove the inductive step. 
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Suppose that for some xq & S, S \ {xq} = ZYi U where Ui and U2 are disjoint and 
Ui is a finite union of finite trees VVj , j = 1, 2, . . . , A;. Assume that Yq = xq and Xq e Wi 
(the argument is analogous if the roles of X and Y are interchanged). Suppose without 
loss of generality that Xq e Wi. Let xi 7^ xq be the vertex of Wi closest to a;o, and 
VVi \ {a^i} = VV2 U W3, where W2 and W3 are disjoint, and W3 is the edge joining xq 
and xi. We will first assume that W2 7^ 0. We will show that for some P3 > and some 
stopping time Tg < 00, with probability greater than ps, we either have d(Xx'g, Ytq) < e 
or Ytq = xi and Xtq E \V2- 

If d(Xo, Iq) < then we can take Tq = 0. 

Assume that d(Xo, Iq) > £• Let X2 G W3 be the point with d(xo, 2:2) = e/3. Note that 
d(a;2, lo) < (l/2)d(x2, Xq). Let T7 = inf{t > : = a;2} and Tg = inf{t > : = X2}. 
By Lemma 3.4, for some s > 0, 

P(T8 < > P(T7 < s). 

Hence, with probability ^3 > 0, Tg < T7 and either X and Y have met by the time Tg, 
or X is on the opposite side of Irg in <S than xq. Let Tg = inf{t > : Yt = xi] and 
Tio = inf{t >T^:Yt—xo}. Since «S contains only a finite number of finite trees 7^, there 
is an upper bound on the edge length in any tree 7^, say, p < 00. This and the fact that 
d(x2,a;o) = s/3 imply that P{Tq < Tiq) > p4 for some ^4 > that may depend on e. If 
the events {Tg < Tj} and {Tg < Tiq} hold then either X and Y have met by the time Tg 
or Yxg) = xi and Xt^) G W2. 

We note that if VV2 = (i.e., Wi is a single edge) then the same argument proves 
that for some ^3 > and some stopping time Tq < 00, we have d(Xj'g,lT6) ^ ^ with 
probability greater than psp^. 

(iii) Let us rephrase the claim proved in step (ii). We have shown that for some 
P5 = P3P4: > and some stopping time Tg < 00, with probability greater than ps, we 
either have d(Xj'g, l^e) < £ or Yt^ and Xxg satisfy the same assumptions as Yq and Xq, 
but relative to graphs Ui = W2 and V(2 = S \ ({^2} U W2) in place of Ui and W2. Recall 
the claim proved in part (i) of the proof and the final remark in step (ii). Note that I4i 
has at least one edge less than Ui so by induction, we can repeat the inductive step (ii) a 
finite number of times and show that with a probability pe > 0, X and Y come within s 
of each other before some time ti < 00. It is easy to check that pq and ti can be chosen 
so that they do not depend on the starting points of X and Y. The Markov property and 
induction can be used to show that X and Y have to come within e of each other by the 
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time jti with probability greater than 1 — (1 — pe)-' . We let j ^ oo to see that X and Y 
come within e of each other at some finite time with probability one. Since e G (0,ro/3) 
is arbitrary, the coupling is not shy. □ 



Example 3.6. Suppose that 5 is a tree with the property that it has a "backbone" that 
is topologically a line, with a finite or countable number of finite trees attached to it. See 
Fig. 3.4 for an example. Recall that we have assumed that each edge has length at least 
ro > 0. 





Figure 3.4. 



Let . . . , x-2, X-i, xq, xi, X2, ■ ■ ■ be the sequence of points along the "backbone" U 
where the side trees are attached (the sequence can be finite or it can extend to infinity in 
one or two directions). If the sequence extends to infinity in both directions and the graph 
is invariant under a non-constant shift, then, according to Example 3.3, there exists a shy 
coupling. 

Assume that 

(i) the diameters of the side trees are uniformly bounded and 

(ii) {xk} does not extend to infinity in both directions, or {xk} extends to infinity in both 
directions but the family {xk} is not shift-invariant, i.e., for every c 7^ 0, there exists 
Xj such that xj + c Xk for all k. 

We will show that under these assumptions there is no shy coupling. 

Let Q be the "projection" of S on U, i.e., Q{x) = x for x E and Q{x) = Xfc if a; 
belongs to a tree which is attached to U at Xk- We will identify the "backbone" U with 
the real line so that we can think of Q{X) and Q{Y) as real-valued processes. Let X and 
y be a coupling of Brownian motions on S. If X makes an excursion into a side tree 
then Q{X) remains constant on the excursion interval (including the endpoints). Hence, 
Q{X) and Q{Y) are continuous processes. It is easy to see that they are local martingales. 
Informally speaking, they are Brownian motions frozen on some random intervals. The 
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process Zi = Q{Xt) — QiYt) is also a continuous local martingale. Suppose without loss 
of generality that > and let 



To = inf{t > : = 0}. 



Then Zt^m is a non-negative local martingale and so it must have an almost surely finite 
limit Zoo on {Tq = oo}. 

We will show first that P(To = oo) = P(7o = oo and Z^o = 0), in other words, 
P(To = oo and Zoo > 0) = 0. Let S be a Brownian motion on M with Bq — 0. Consider a 
small £ > and to > 0, let S e (0, e) be such that 



Let Ti be the first time when all of the following conditions hold: d{XTi,V() > e, the 
distance from Xti to any vertex of S is greater than e, Yt^ G W, and infjt d(lTi, Xk) > 6 
(the argument is analogous if the roles of X and Y are interchanged). If Ti < oo then with 
probability pi/2 or greater, X will stay on the same side tree over the interval [Ti, Ti +^0], 
while Y will move away from Yt-^ by more than d units over the same time interval. Hence 
with probability pi/2 or greater, Zt will have an oscillation of size at least 6 over the interval 
[Ti, Ti+to]. We proceed by induction. If < 00 then we define Tk+i — Tio9rp^_^_f^+Tk+tQ, 
where 6. is the usual Markovian shift operator. Then with probability greater than pi/2, 
Zi has an oscillation of size at least 5 over the interval [Tfc+i,Tfc_|_i + to], independent of 
whether that happened over any interval [Tj,Tj +to], J < ^- Hence, with probability one, 
either = 00 for some k or Zt has an infinite number of oscillations of size 8 over disjoint 
intervals of length to, and, therefore in the latter case, Zt does not have a limit as t — > 00. 
Applying the above argument to a decreasing sequence of n > 1} and a decreasing 
sequence of {5^, n > 1} both tending to zero and after deleting a null set from f2, we may 
and do assume that for every a; e and for every £„, there is some N > 1 such that for 
every j > N, with and Sj in place of s and S above, either Tj.{uj) = 00 for some k or 
Zt{oj) does not have a limit as t — > 00. The processes X and Y are recurrent because the 
one-dimensional Brownian motion is. Hence, for every Xj, each one of them will enter the 
side tree attached to U at xj infinitely often. After deleting a null set from Q, we may and 
do assume that the aforementioned property holds for every co & Q. 

For a; e {Tq = 00 and Z^o > 0}, let 




where pi :—P \ sup \Bt\ < e 




c(a;) = lim Zt{u!) > 0. 
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We choose an Xj, relative to c{uj), as follows. If {xn} does not extend to — oo (oo) then 
we let Xj be the leftmost (rightmost, resp.) point of the sequence. Otherwise we fix an Xj 
with the property that Xj + c{u>) ^ Xk for all k (such an xj exists by assumption). Note 
that both X{ijj) and Y{(jj) enter the side tree attached to U at xj infinitely often. Hence 
one can find some e > from {en,n > 1} and S > from {Sn,n > 1}, and an increasing 
sequence of random times {Sk, A; > 1} with limfc^oo Sk = oo such that all of the following 
hold. One has d{Xs^,h() > e, the distance from Xs^ to any vertex of S is greater than 
£5 £ ^) ^iid inf„ d(lsj^, Xn) > 5 for every A; > 1 (or the statement will hold with the 
roles of X and Y interchanged). 

Hence, all stopping times {T}^, k > 1} defined in the proceeding paragraph are finite. 
We have shown that this event implies that Zoo{uj) does not exist. This contradiction proves 
that P(To = 00 and Zoo > 0) = and therefore P(To = 00) = P(To = 00 and = 0). 

Recall that we have assumed that all the edges have length at least ro > 0. So on 
{To = 00 and Zqo = 0}, by the recurrence of the one-dimensional Brownian motion, we 
have liminf^^oo d(Xt, Yf) = 0. We now only need to exam uj e {Tq < 00} and to prove 
X{lj) and Y{uj) will come arbitrarily close to each other. 

Consider any e G (0, ro/4), where ro > is a lower bound for the length of any edge in 
S. We want to show that with probability one, there exists t such that d{Xt, Yt) < e. We 
have already proved that liminft^oo d{XtjYt) = on {Tq = 00}. On {Tq < 00}, at time 
To, either Xtq — Ytq or one of processes {^To7 ^Tq} is at some x^ and the other process 
is in a side tree T attached to U ai Xk- Without loss of generality, assume that Xtq = Xk 
and Yto G If d(XTo, Yt^) < e then we are done. Suppose that d(XTo, iro) > £• Let zq 
be the point at the edge e of T that is attached to U with d{zo,Xk) = d{zo,U) = e/4. Let 
be the first time after Tq when Xt = zq. By Lemma 3.4, with probability p2 > 0, X^ 
reaches zq after Tq before Yt gets there. If this event occurs, both and Ys^ will have 
distance at least e/4 away from U. Let Ri be the first time after Si when both Xt and Yt 
are outside T. An argument analogous to that in parts (ii) and (iii) of Example 3.5 shows 
that with probability ps > 0, the processes X and Y will meet during the time interval 
[Si, Ri]. In other words, conditioning on {Tq < 00} and d(XTo,lTo) > with probability 

df 

at least = P2P3 > 0, the processes X and Y will meet between times Tq and the first 
time Ri when they are both outside T. We define for A; > 2, 



Sk = Sio dR^_^ + Rk-i and Rk = Ri o 63^ + Sk- 
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By the strong Markov property of {X, Y) , 

p( inf diXt,Yt)>e] <(1-P4)'. 



Letting ^ oo, we get 



P inf d{Xt, Yt)>e] <0 

\te[o,oo) ' 



for every e > and thus {X, Y) is not a shy couphng. 



Example 3.7. Suppose that S is composed of a loop U with a finite number of finite trees 
attached to it at points Xk, and the family {xk} is not rotation invariant in the following 
sense. We can assume without loss of generality that U is isometric to the unit circle. For 
every c 7^ 0, there exists Xj such that xje^^ 7^ Xk for all k. See Fig. 3.5 for an example. We 
will show that in this case there is no shy coupling. 




Figure 3.5. 



Our argument will be very similar to that in Example 3.6. Recall the "projection" Q 
from the previous example. We have Q{x) = x for x eU and Q{x) = Xk if x belongs to a 
tree that is attached to U at x^. Hence, Q{Xt) may be regarded as a continuous process 
on the unit circle. We now choose a (random) continuous function Ox '■ [0, 00) ^ M so 
that Q{Xt) = e*®^*-*-* for all t > 0, in the complex notation. We define 0y in an analogous 
way. Note that Gx and Qy are martingales. Therefore, Zt = Ox(^) — 6)y(^) is also a 
martingale. We can now repeat the argument from Example 3.6 to show that there does 
not exist a shy coupling. □ 



Example 3.8. Examples 3.3 and 3.5-3.7 may appear to suggest that if a graph has a vertex 
with degree 1 then a shy coupling exists only if there exists an isometry of I : S ^ S with 
no fixed points. We will show that this is not the case. Our example is illustrated in 
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Fig. 3.6. In this case, every isometry I : S ^ S has a fixed point. Nevertheless, we wiU 
show that there is a shy couphng in S. 




We will describe below the transition mechanism for {X, Y) on some random intervals 
of time. We will assume that the transition probabilities of {Y,X) are the same as those 
of (X, y). Hence, there is no need to describe cases symmetric to those discussed below, 
in the sense that the initial positions of X and Y are interchanged. 

Suppose that all edges Ai, A2, . . . , A7 have the same length, say 1. We assume that 
Xq = X2 and Yq = xs- 

(i) Suppose that for some stopping time Ti we have Xt^ = X2 and = x^. Then 
we let T2 = inf{t >Ti:Xt^ {Ai U A2 U A3) \ {xi}}. Let / : (Ai U A2 U A3) \ {xi, X3} ^ 
(A3 U A4 U A5) \ {x2, X4, X5} be the one-to-one isometry satisfying /(X2) = X3, I{Ai) = A3, 
/(A2) = A4 and /(A3) = A5. We let Yt = I{Xt) for t e [Ti,T2]. Note that at the stopping 
time T2, we have one of the following configurations of the two particles: {Xt2,Yt2) — 
{xi,X4), or (Xt^.Yt^) = (xi,X2), or (Xt^.Yt^) = (x3,X5). 

(ii) Suppose that for some stopping time T3 we have Xt^ = xi and Yt^ = X4. Let 
{Xf, t G [T3, T4]} be Brownian motion on S independent of the past with Xto, = a^s, where 
T4 = inf{t > T3 : = X2}. For t G [T3,T4], we let Yt G A4, with diYt.x^) = d{Xt,xi). 
Note that Xt^ = X2 and = X3. 

(iii) Suppose that for some stopping time T5 we have Xt^ = xi and Yt^ = X2- Let 
{Yt, t G [T5, Te]} be Brownian motion on S independent of the past with Yt^ = X2, where 
Tg = inf{t > T5 : 1^ = xi or X3}. We label excursions of Y from X2 that stay in ^.3 
with marks "1" or "2", with equal probabilities, in such a way that the label of any 
excursion is independent of all other labels. Then we let Xt be defined for t G [T5,T6] by 
d(Xt,xi) = d(yt,X2) and the following conditions. If Yt G Ai then Xt G A2, if Yt G A2 
then Xt G Ai, if Yt G A3 and t belongs to an excursion marked "1" then Xt G Ai, and 
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if Yt e ^3 and t belongs to an excursion marked "2" then Xt e A2. At time Tq we have 

(XTe,YTe) = ix2,Xi) or (Xts^Yt^) = {X2,X3). 

Note that S is symmetric with respect to the hne containing A4. If for some stopping 
time T7 we have (Xt^, IT7) = {x3,X5), or (Xy^, IT7) = {^4, xq), or (Xy^, IT7) = {^5, xe), or 
one of these conditions is satisfied with the roles of X and Y interchanged, then we define 
the coupling on an appropriate random interval in a way analogous to that in (i)-(iii), 
using the symmetry of S. 

The above definitions for the "local" behavior of the coupling and the strong Markov 
property can now be used to define a process (X^, Yf) for alH > 0. It is easy to see that 
the stopping times analogous to Ti , T3 and T5 will not have a finite point of accumulation. 
It is also easy to check that almost surely d(Xt, Yt) — 1 for every t > 0. □ 

4. Reflected Brownian motion in Euclidean domains. 

This section is the closest in spirit to the papers and problems which inspired the 
present research project. Suppose that D C M.^ is a bounded connected open set which 
is either convex or has a boundary. We will consider couplings (X, y) of refiected 
Brownian motions in D, defined as follows. Let n{x) denote the unit inward normal vector 
at a; e dD. Let B and W be standard planar Brownian motions with Bq = Wq = defined 
on the same probability space and consider the following Skorohod equations, 



Here L is the local time of X on dD, i.e., a non-decreasing continuous process which 
does not increase when X is in D: l]j{Xt)dLf = 0, a.s. Equation (4.1) has a unique 
pathwise solution (X, L^) such that Xt E D for all t > (see [Ta] when D is convex 
domain and [LS] when D is ). The "reflected Brownian motion" X is a strong Markov 
process. We point out that B is uniquely determined by X, and vice versa. The same 
remarks apply to (4.2), so, as a pair, (X, Y) is also strong Markov. 

For a continuous semimartingale M, the symbol (Z) will stand for its quadratic vari- 
ation process. When M = (M\ • • • , M'^) and AT = (Z^, • • ■ , Z'^) are two continuous R'^- 
valued semimartingales we will use {M,N) to denote Ylt j=i{^^i^^)- Note that the 
matrix-valued process {{M'^, N^))i<ij<d is non-negative definite and so t 1— > {M,M)t is 
always non-decreasing. For a,b e R'', we use a ■ b to denote the inner product between a 




(4.1) 



(4.2) 
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and b. We will use d{x,y) and \x — y\ interchangeably for the Euclidean distance between 
x,ye R'^. 



Theorem 4.1. Assume that D cW^ is a bounded convex domain. Let X and Y be two 
reQecting Brownian motion on D given by (4.1)- (4. 2). 
(i) Suppose that there is a strictly increasing function (p with (p{0) = such that 

d{\X - > ^{\Xt - Yt\)dt for t < ao, 

where ctq := inf {t > : = Y^}. Then (X, Y) is not a shy coupling, 
(a) Suppose D is strictly convex. Assume that {X — Y,X — Y)t (this is the same as 
{B — W,B — W)t) ) has a sublinear growth rate as t ^ oo, that is, 

lim {X -Y, X -Y)t/t = almost surely. 

Then (X, Y) is not a shy coupling. 

Proof. Note that 

Xt-Yt = Xo-Yo + {B-W)+ [ n(X,)dLf - / n{Y,)dLl 

Jo Jo 

is a semimartingale. Define Rt := |X^ — l^tp. By Ito's formula, 

dRt = 2(Xt - Yt) ■ diXt - Yt) + d{X-Y,X- Y)t 

= 2iXt - Yt) ■ d{Bt - Wt) - 2{Yt - Xt) ■ n{Xt)dLf - 2(X, - Yt) ■ n{Yt)dLY 
+ d{X -Y,X -Y)t. (4.3) 

(i) Let a > be a constant whose value will be chosen in a moment and f{r) — — r~" 
for r > 0. Then f'{r) = ar~"~^ > and f"{r) = {—a — l)ar~"~^ < for r > 0. Define 
Ut := f{Rt) = f{\Xt - Yt\^). By Ito's formula, we have 

dUt = f'{Rt)dRt + ^-f"{Rt)d{R)t = dMt + dVt, 

where 

dMt = 2aRt''-^2{Xt - Yt) ■ d{Bt - Wt) 

and 

dVt = -2aRt''-' {{Yt - Xt) ■ n{Xt)dL^ + {Xt - Yt) ■ n{Yt)dLl) 

+ aRt''~'^d{X -Y,X-Y)t- 2a{a + l)Rt''~'^d{\X - Yf)t (4.4) 
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are the local martingale and bounded variation parts, respectively. We claim that for every 
£ > 0, 

Te :=inf{t>0: \Xt-Yt\<e) 

is finite almost sure. Suppose that P(Te = oo) > for some £ > 0. We will show that this 
leads to a contradiction. 

Since D is a convex domain, for e dD, n[x) is well defined and 

{y — x) • n(a;) > for every y & D. 

Note that the local time (respectively, L^) does not increase when X (respectively Y) 
is on a subset of dD having zero Lebesgue surface measure. Hence, (4.4) yields 

dVt < aRt^-^diX -Y,X-Y)t- 2a{a + l)Rt''-^d{\X - F p)^. (4.5) 

Note that d{X -Y,X-Y)t = d{B -W,B- W)t < Mt. This, (4.5) and the hypothesis in 
part (i) of this theorem imply that on {Tg = oo}, 

dVt < AaRt^-^dt - 2a{a + l)Rt''-^(p(s)dt 

< -2aR-''-'^ ((a + !)(/?(£) - 2Rt) dt 

< -2a£-''-^ ((a + l)(p{e) - 2diam(D)) dt. 

For a fixed £ > 0, we can find a > sufficiently large so that for some A > 0, 

dVt < -Xdt for every t > on {T^ = oo}. (4.6) 

The continuous local martingale M is a time change of Brownian motion. By the law 
of iterated logarithm for Brownian sample path, for almost all u) G {Tg = cxd}, there is an 
unbounded increasing sequence {tfc, k > 1} such that sup;j>i \Mti^{u>)\ < oo. This and (4.6) 
imply that Ut^{Lv) = Vo{u) + Mt^{u)) + Vti^{co) tends to — oo as ^ oo on {Tg = oo} a.s. 
Consequently, \Xt^ — Y-tJ goes to as /c — >^ oo on {Tg = oo} a.s., which is a contradiction. 
This proves that particles X and Y come arbitrarily close to each other in finite time and, 
therefore, X and Y is not a shy coupling. 

(ii) Now assume that D is bounded and strictly convex and {X — Y,X — Y)t has a 
sublinear growth as t — > oo. The strict convexity implies (in fact, it is equivalent to) the 
following condition. For every small £ > 0, there is a constant > such that 

(y — x) • n{x) > as\x — y\ for every x e dD and y & D with |a; — y| > s. (4.7) 
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Let a denote the surface measure on dD. Since reflecting Brownian motion in D is a 
recurrent FeUer process, it foUows from the Ergodic Theorem that 

<j{dD) , 

hm = = lim almost surely. 

t— >oo t 2 LD t— >oo t 

For every £ > 0, define := mf{t > : \Xt — Yt\ < e}. On {Tg = oo}, we have from 
above and (4.7) that 

limmf 1 (^-2 ^ {Yt - X^) ■ n{Xt)dLf - 2 ^ (X, - F) • n{Yt)dLj + {X-Y,X- Y)^ 

< liminf - (~2£a,Lf - 2eaeLj + {X-Y,X- Y)t) = _?£^£^iM < q. (4.8) 
t^oo t \D\ 

On the other hand, := 2 jl{Xt — Yt) ■ d{Bt — Wt) is a continuous martingale and thus 
is a time-change of one-dimensional Brownian motion. By the law of iterated logarithm 
for Brownian sample path, for almost all uj e {Tg = oo}, there is an unbounded increasing 
sequence {tk,k > 1} such that sup^>i \Mt^^{uj)\ < oo. This, (4.3) and (4.8) imply that 
lim^^oo Rt = — a-S- on {T^ = oo}. Since Rt > 0, we conclude that P{Ts = oo) = for 
every e > and so (X, Y) is not a shy coupling. □ 



In the remainder of this section, we take d — 2, but this is only for notational con- 
venience. We will show in the next example that the method of proof of Theorem 4.1, 
based on the Ito formula, does not extend to arbitrary couplings. The example may have 
some interest of its own. We will show in Theorem 4.3 below that, in fact, there is no 
shy coupling of reflecting Brownian motions on any bounded C"'^-smooth strictly convex 
domain. 

Example 4.2. We will show that there exist planar Brownian motions B and W with 
the property that d{Bt, Wt) = ^J2t + d(So, W^Y for t > 0, assuming that Bq ^ Wq. In 
particular, the distance between the two processes grows in a deterministic way. 

Suppose that {B, W) has the above mentioned property with Bq and Wq taking values 
in D. Let X and Y be the pathwise solutions of (4.1)-(4.2) but with the above B and W 
in place oi xq + B and yo + W there. We have 

d{\Xt-Yt\^) = d{\Bt-Wt\^) = 0, 

while 

d{X -Y,X-Y)t = d{B -W,B-W)t = d {\Bt - Wt]'') = 2t. 
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So Vt in (4.4) becomes 



dVt = -2ai?,-""' {{Yt - Xt) ■ n{Xt)dLf + {Xt - Yt) ■ n{Yt)dLj) 
+ 2aRt''-^dt. 

Hence the method used in the proof of Theorem 4.1 does not work for this couphng (X, Y). 
Moreover, since jX^ — y^] grows deterministically when both Xt and Yt are away from the 
boundary and decreases when one of them is on the boundary, neither \Xt —Yt\ nor any 
deterministic monotone function of |Xt — Yi| is a submartingale or a supermartingale. 

We now present the construction of B and W with the properties mentioned above. 
Let 5 be a Brownian motion in starting from a^o- For a vector v = (a, 6) e IR^, we 
use v-^ to denote its orthogonal vector (6, —a). Let yo G be a point different from xq. 
Consider the following SDE for W in with Wq = yo: 

dWt = 1^.,, ^ ^ „ {{{Wt - Bt) ■ dBt) {Wt - Bt) - {{Wt - Bt)^ ■ dBt) {Wt - Bt)^) . 

In words, at any given time t > 0, takes a synchronous step with Bt along the direction 
Wt—Bt, while Wt moves in the opposite direction but with the same magnitude as Bt along 
the perpendicular direction (Wt — Bt)-^. The above SDE for W has a unique solution up to 
T := inf{t > :Wt = Bt}, since the diffusion coefficients are C°° up to that time. It can be 
computed directly that d (jWt — -B^p) = 2dt and consequently |Wt — = |a;o — + 2t. 
So T = oo. It is standard to check that W = (VF^, W"^) is a continuous local martingale 
with {W\ W'^)t = t for z = 1, 2 and {W^, W'^) — 0. Therefore W \s b, Brownian motion in 
starting from yo- □ 

We will show next that in a C-'^-smooth strictly convex domain D, every coupling of 
reflecting Brownian motions on D must come arbitrarily close to each other in finite time. 

Theorem 4.3. Suppose that D is a bounded convex planar domain with a -smooth 
boundary that does not contain any hne segments. Then there does not exist a shy couphng 
{X, Y) of reEected Brownian motions in D. 

Proof. The idea of the proof is inspired by differential games of pursuit (see [F] ) . We will 
show that with positive probability, one of the particles will pursue the other one in such 
a way that the distance between the two particles decreases either because the diffusion 
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component of the second process does not move the second particle sufficiently fast or the 
second particle hits the boundary and is pushed back towards the first one. 

Step 1. We will define several constants Sk in this step. The definitions will be labeled 
(a), (b), (c), etc. Each of these definitions is really a simple lemma asserting the existence 
of a constant with stated properties. Since the proofs do not need more than high school 
geometry, we omit most of the proofs. The constants Sk are defined relative to each other, 
but £k may depend only on the values of Sj for j < k. 

(a) Let £o > be so small that for every x & D with d(a;, dD) < eq there exists a 
unique point in dD whose distance from x is minimal. 

We make £o > smaller, if necessary, so that the following is true. Consider any point 
y G dD and let CSi be the orthonormal coordinate system such that y = G dD and 
n(0) lies on the second axis. Write n(a;) = (ni(a;), n2(a;)). Then |ni(a;)| < n2(a;)/100 for 
X e dDnB{0,eo) in CSi. 

We fix an arbitrary £i e (0, £o]- It will suffice to prove that for any xo,yo e D, if 
(Xo, Yo) = {xq, yo) then, with probability one, there exists t < oo such that d{Xt, Yf) < e\. 

(b) The angle between two vectors will be denoted Z(- , •), with the convention that 
it takes values in (— 7r,7r]. Since D is a bounded and strictly convex domain, there exists 
£2 e (0,7r/2) such that for every x G dD and y & D satisfying d{x^y) > £i/2. 



(c) Let r) be the cone spanned by {n{y), y G dD (1 B{x,r)}. Since dD is C^- 
smooth, £(x,r) is a wedge. Hence, all linear combinations of vectors in £{x,r) with 
non-negative coefficients belong to C{x,r). An easy approximation argument shows that 
if Xt G B{x, r) for aU t G (s, u) then n{Xt)dLf G C{x, r). 

We will now choose £3 G (0, £1/8). Consider vectors v and w satisfying the following 
conditions, relative to Xcyo? and £3. 

If d(a:o, dD) < £3 then v G C{xq, 26^) and |v| < £3. If d(xo, dD) > £3 then v = 0. 

If d(yo, dD) < £3 then w G £(yo, 2£3) and |w| < £3. If d(yo, dD) > £3 then w = 0. 

We will show that (4.9) implies that we can find sufficiently small £3 > so that 
the following is true. Suppose that XQ,yo G D with d{xo,yo) > £1. Assume that xi G 
;B(a;o,2£3) and yi G B{yo,2e3). Then 



l{n{x),y-x) G [-7r/2 + £2, 7r/2-£2]. 



(4.9) 



d(a;i + V, yi + w) < d(a;i, yi). 



(4.10) 



To see this, choose some X2 and y2 so that the following conditions hold. 
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If d{xo, dD) < £3 then X2 £ dD fl B{xq^ 26^) with v = cr\{x2)- If d(a;o, dD) > £3 then 

X2 = Xq. 

If d{yo, dD) < £3 then y2 G 5-D fl B{yo, 2£3) with w = cn(j/2)- If d(yo5 ^-D) > £3 then 

y2 = yo- 

Since £3 < £1/8, 

|a^2 - y2\ > \xo - yo\ - |a;o - a;2| - \yo - 2/2! > £1 - 4£3 > si/2. 
Thus by (4.9), we have for sufficiently small £3 > 0, 
d(a:i + V, 2/1 + w)^ 

< \xi - + |v - w|^ + 2{xi - yi) • (v - w) 

< \xi - yiP + |v - wp + 2{x2 - 2/2) • (v - w) + 2(1x1 - X2\ + \yi - y2|)|v - w| 

< \xi - + |v - - 2(sin£2)|a;2 - ?/2|(|v| + |w|) + 4(8£3)|v - w| 

< \xi - + (|v| + |w|) (6£3 - 2(sin£2)(£i/2) + 32£3) 

= \xi - yip - (|v| + |w|) (£1 sin£2 - 38£3) (4.11) 

< d{xi,yif. 

(d) Since D is bounded, we can find £4 > and N < 00, such that if xi,X2, ... is a 
sequence of points with xi e D, d{xk,Xk-i) > £3/8 and \/-{xf~ — Xk-i,Xk+i — Xk)\ < 64^ 
for all k then xn ^ D. 

(e) We choose £5, £6 > so that the following is true. Suppose that xo,yo,xi,yi and 
X2 satisfy the conditions d(a;o,yo) ^ ^i; d.{xi,yi) > £1 and 

\l{xi-yi,xo-yo)\ < £5, 

d(a;i, xo + (£3/4)(?/o - xo)/d{xo, yo)) < £5, 

d(a;2,a;i + (£3/4)(yi - xi)/d{xi,yi)) < ee- 

Then \/.{xi — xq, X2 — xi)\ < £4. 

We make Sq smaller, if necessary, so that Eq < £3/8. 

(f) We make £6 > smaller, if necessary so that the following holds. Suppose that 
xo,yo,xi and yi satisfy the conditions d{xo,yo) > Si, 

d(a;i,a;o + (£3/4)(?/o - a;o)/d(a;o, yo)) < ^6, (4.12) 
d(yi,yo + (£3/4)(yo - a;o)/d(a;o,yo)) < se- (4.13) 
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Then \l{xi -yi,xo-yo)\ < £5- 

(g) We can find £7, Sg > with the following properties. Suppose that xq, yo e D, 
xi,X2 e IR^, d{xo,yo) > £1 and the following conditions are satisfied. 

If d(a;o, dD) < £3 then v e jC{xo, 2£3) and | v| < £3. If d{xo, dD) > £3 then v = 0. 

If d{yo, dD) < £3 then w e jC{yo, 2es) and |w| < £3. If d{yo, dD) > £3 then w = 0. 

Assume that 



Then d(a;i + v, yi + w) < d(xo, yo) - e-j- 

(h) It is easy to sec from (4.11) that we can strengthen (4.10) as follows. We can 
make £7 > smaller, if necessary, so that if |v| > £6/2 or |w| > £5/2, and the assumptions 
stated in Step 1(c) hold then 



Ste]) 2. Suppose that Xq, Yq ^ D with d(Xo, Yq) > ^i- Consider the following events. 



Fi(t) = {d{Xt,Yt) < d(Xo,yo) -£7}, 
F^it) = {\liXt - Yt, Xo -Yo)\< £5}, 

Fsit) = {d{Xu Xo + {e3/4){Yo - Xo)/d{Xo, Yo)) < se A ss}, 
F^{t) = {d{Xt,Yt) < d{Xo,Yo)+er/{m)}, 
F^it) = (Flit) U F^it)) n Fsit) n F4(t). 



We will show in Step 4 that P(F5(ti)) > pi for some > that do not depend on Xq 

and Yq. 

Let £9 = Sr/{16N) A £3/8 A £6/8 A £3/5. Recall that B and W are Brownian motions 
with Bq = Wq = driving X and Y in the sense of (4.1)- (4.2) and let 



d(a;i, xo + (£3/4)(|/o - xo)/d{xo, yo)) < Sg, 
d{yi,yo + {es/A){yo - xo)/d{xo,yo)) > ee/2, 
d(yi,2/o) < £3/4 + £8- 



(4.14) 
(4.15) 
(4.16) 



d{xi + V, yi + w) < d{xi, yi) - le-j. 



(4.17) 



^i(0 



{Bt e B ((£3/4)(Fo - Xo)/d{Xo, Yo), £9)} , 




sup \Bs - {s/t)Bt\ < £9 

s6[0,t] 



{|m| <£3/4 + £9}, 




<^ sup \Ws-{s/t)Wt\<eQ) 

[s€[0,t] J 

Ai{t)nA2{t)nAs{t)nA^{t). 



A^{t) 
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We will argue in the rest of this step of the proof that P(^5(ti)) > pi for some pi,ti > 0. 
In later steps, we will show that ^5(t) C -F5(t). 

Recall that 5 is a two-dimensional Brownian motion with Bq = and let = inf{t > 
: \Bt\ > r}. Note that, by Brownian scaling, 

P(Tr <t) = P I max \Bs\ > r | = P ( max IBJ > r/Vi] . 

\0<s<t ) \^<s<X ) 

By the large deviations principle (see [RY], Ch. VIII, Thm. 2.11), 

lim (2t/r^)logP(T^ < t I So = 0) = -1. (4.18) 

Let 

ro = £3/4, 
A^it) = {Tr, < 0, 

A7 = {Bt^, e Biies/4)iYo-Xo)/d{Xo,Yo), £9/2)}, 

A8{t) = \ sup \Bs-BTj<eg/2 

Clearly ^6(^), A7 and Agit) are independent, and ^6(*) n ^7 n Agit) C Ai{t). So 

P{A,{t)) > PeiA^it)) ■ P{Ar) ■ PiAsit)). 

Note that limt^o P{^8{i)) = 1 by the strong Markov property applied at T^^ and the event 
At is independent of t. This together with (4.18) yields 

liminf tlogP(Ai(t)) > -rl/2 = -(£3/4)2/2. 

We obtain directly from (4.18) that 

limsup tlogP(A^(t)) < -(£3/4 + £9)^2. 

This implies that for all sufficiently small t > we have 

P(^i(t)n^3(*))>0. (4.19) 

The process {Bg — {s/t)Bt, s e [0,t]} is Brownian bridge with duration t seconds, i.e., 
Brownian motion starting from and conditioned to be at at time t. If Tr denotes the 
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hitting time of r by the absolute value of the Brownian bridge then we have a formula 
analogous to (4.18), limsupi/^2^o(2*/^^) logP(7; < t) < -1. Thus 

limsupnogP(^^(t)) < -£9/2. (4.20) 

Note that Brownian motion B = {Bt,t > 0} is a Gaussian process. Since for every 
< s < t, Bg — {s/t)Bt and Bt have zero covariance, the process {Bg — {s/t)Bt, s e [0, t]} 
is independent of Bf. Similar remarks apply to W. Hence 

F{A,{t) n A2{t)) = P(Ai(t))P(A2(t)) and F{A^{t) n A^{t)) = P{As{t))P{A^{t)). 

This, (4.19) and (4.20) imply that 

P{A5{ti)) > pi for some h.pi > 0. (4.21) 

Step 3. Let Rt = /o n{Xs)dLf . We will show that if (to) 0^2(^0) holds for some to > 
then \Rs\ < Asg for every s G [0,to]- 

Since Ai{to) and ^2(^0) hold, \Bt\ < £3/4 + 2eg < £3/2 for every t < to. Thus when 
d(Xo, dD) > £3, Xo + Bg E D for all s e [0, to]- By the uniqueness of the solution to (4.1), 
Xg = Xo + Bg for s e [0, to] , and Rg = for all s e [0, to] • 

Suppose that d{Xo,dD) < £3. By the assumptions made in Step 1(a), there exists 
a unique point y e dD with the smallest distance to Xq. Let CSi be the orthonormal 
coordinate system such that y = E dD and n(0) lies on the second axis. Recall from Step 
1(a) that n{x) = (ni(a;), n2(a;)) and |ni(a;)| < n2(a;)/100 for x e dDnB{0,3e3) in CSi. 
Write Rt = {Rl,Rt)- By the opening remarks in Step 1(c), Rt E C{0,r) if Xg e B{0,r) 
for all s<t. This implies that \Rl\ < i?f/100 if Xg e B(0,3£3) for all s<t. 

Let Ti = inf{t > : |-Rt| > 4£9}. We will assume that Ti < to and show that this leads 
to a contradiction. Since \Bt\ < £3/2 for every t < to, we have \Bt + Rt\ < £3/2 + 4£9 < £3 
for every t < Ti. Hence, 

\Xt\ < \Xo\ + \Bt + Rt\ < £3 + £3 = 2£3 for every t < Ti. 

It follows that Rt > for every t < Ti and \R^^\ < R'^JlOO. Note that, by Step 
1(a), the slope of the tangent line at points in B(0,3£3) fl dD is between —1/100 and 
1/100 and that d{XTi,Xo + Bt^) = \Rti \ = 4£9. The last observation and the fact that 
e dD imply that d(Xo + BT^,dD) > (2/3)|i?Til = 8£9/3. Since ^i(to) and ^2(^0) 
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hold, d(Xo + Bt,D) < 2eQ for every t < to and in particular for t = Ti. This contradiction 
proves the claim that |-Rs| < 4£g for every s e [0,to]- 

Let Rt = n{Ys)dLY . We claim that if ^3(^0) H ^4(^0) holds for some to > then 
I-Rsi < £3 for every s G [0,to]- To see this, observe that d{Yo + Wt,D) < £3/4 + 269 < £3/2 
for t <to and use that same argument as in the case of Rf. 

Step 4- Fix to > 0. We will show that ^5(^0) C -^5(^0)- Assume that ^5 (to) holds. 
Since Ai{to) holds, we have in view of Step 3, 

diXt„Xo + (£3/4)(yo - Xo)/d(Xo,yo)) 

< d(S,„, (£3/4) (yo - Xo)/d{Xo, Yo)) + \Rt,\ < £9 + 4£9 < A Ss- 

In other words, ^3(^0) holds. 

Since Ai{to) and A3 (to) hold we have, using simple geometry, 

d(Xo + Bt„Yo + Wt,) < d{Xo, Yo) + 3£9. (4.22) 

We will apply (4.10) with xi = Xq + Bt„, yi = Yq + VFt;,, v = Rt^ and w = Rt^. We have 
assumed that ^5(^0) holds so \Bt\ < £3 and < £3 for t < to- This and Step 3 imply 
that for t <to, 

d{Xt,Xo) < \Bt\ + \Rt\<2e3, 
and similarly d{Yt,Yo) < 2ez. By the opening remarks in Step 1(c), 

V = / ° n{Xt)dLf e CiXo, 283) and w = / ° n{Yt)dLj e C{Yo, 283). 
Jo Jo 

By Step 3, |v| < £3 and |w| < £3. We have shown that all the conditions listed in Step 

1(c) are satisfied so we can apply (4.10) to obtain 

d(x^„,y^j<d(Xo,yo) + 3£9. 

This proves that -F4(to) holds. 

It will now suffice to show that if -F2(to) does not hold then Fi(to) does. Assume that 
^2(^0) does not hold. 

If all of the following conditions hold, 

\Bt, - (£3/4)(yo - Xo)/d(Xo, yo)| < £6/2 A C8, (4.23) 

\Wt, - (£3/4)(yo - Xo)/d{Xo, Yo)\ < £6/2, (4.24) 

d{Xto,Xo + Bto)<ee/2, (4.25) 

d(y^„,yo + mo) <£6/2, (4.26) 
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then (4.12) and (4.13) hold with {xi,yi) = {Xt^,YtJ and {xo,yo) = {Xo,Yo), and this 
imphes F2{to), which is a contradiction. Hence, at least one of the conditions (4.23)-(4.26) 
must fail. The first of these conditions holds because Ai{to) is true. By Step 3, (4.25) 
holds. 

Suppose that (4.26) fails. In view of (4.22), we can apply (4.17) to xi = Xq + Btg, 
yi =Yo + Wto , V = Rto and w = Rt^ to obtain 

diXt^,Yj < d(Xo, Yo) + 3£9 - 287 < d(Xo, Yq) - 87. 

Hence, we have -Fi(to) iii this case. 

Suppose that (4.24) fails. Then, in view of Step 3, (4.14)-(4.16) hold with (xq, j/o) = 
{Xq, Fo), {xi, yi) = {Xq + -Bto, Fo + W^to)' ^ = Rt^ and w = Rt^ and we have 

d{Xt,,Yt,)<d{XQ,YQ)-87. 

Hence, Fi(to) holds. This proves that ^5(^0) C -^5(^0)- 

Step 5. Fix some £1 e (0,£o) and let Sj's be defined relative to £1 as in Step 1. Let p be 
the diameter of D and let A^o be an integer greater than 4p/£7. Recall from (4.21) in Step 
2 that for some > we have P(A5(ti)) > pi. Let 

Si = inf{t > : d{Xt,Yt) < 81 or F5{t) holds} A {2ti). 

By (4.21) and Step 4, P{Si < ti) > pi. Recall that 9 stands for the usual Markov shift 
and define 

^0 = and Sk = Sio es,_, + Sk-i for k>l. 

Recall integer N defined in Step 1. By the strong Markov property, with probability no 
less than pl^^° > 0, we have Sk - Sk-i < h for all k < 2NNq. 

We will argue that if nA;<2iVJVo {"^^ ~ Sk-i < h} holds then d{Xt,Yt) < 8\ for some 
t < 2NNQti. Assume otherwise. Then F5{Si) o e's^^,, holds for every k < 2NNq. In 
particular, F4(^i) o holds for every k < 2NNo. Let Fe{t) = F2{t) n Fz{t). Since 

F^{t) C F^{t) n {Fi{t) U {F2{t) n Fs{i))), for every k < 2NNq, at least one of the events 
Fi{Si) o es,_, and ^6(51) o es,_, holds. 

Consider any j < Nq. If ^6(^1) o es^_, holds for k = 2jN, 2jN + 1, . . . , 2(j + 1)N - 1 
then Xs^. 's and Ys^ 's satisfy the following conditions for k = 2jN, 2jN+l, . . . , 2(j+l) A^— 1, 

d{Xs,^,,Xs, + {8s/4){Ys,-Xs,)/d{Xs,,Ys,))<8e, 
\l{Xs,-Ys„Xs,^,-Ys,^,)\<8,. 
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This implies, by Step 1(e), that for k = 2jN, 2jN + 1, . . . ,2{j + 1)N - 2, 



I ^ {Xsk+i - ^Sk 1 - ^Sk+i ) I < £4 • 

Since ^3(51) o 9sf._^ holds, we also have d{Xs^^^, Xs,^) > £3/8 for the same range of k 
(to see this, recall from Step 1(e) that Eq < £3/8). Hence -^S2(j_,_i);v-2 n^^st be outside D, 
according to the definition of £4 and N in Step 1(d). Since X always stays inside D, at 
least one of the events FQ{Si)o9sf._^ must fail for some 2jN < k < 2{j + 1)N — 1. Hence, at 
least one event Fi{Si)o9s^ _^ holds for some 2jN < k < 2{j + 1)N — 1. Since F4{Si)o9s^_-^ 
holds for every k < 2NNo, there is a reduction of at least £7/2 in the distance between X 
and Y on every interval [S2jN , S2[j+i)N], that is, 

d (^52,w> Ys^jn) < d {Xs,^._,^^, >S2o_i);v) - y for every j e {!,■■■ , Nq}. 
Summing over j we obtain 

d (X5,^„^, < d(Xo, Fo) - ^ < d(Xo, Yo) - 2p < 0. 

This contradiction proves our claim that 

if pl {Sk - Sk-i < h} holds, then d{Xt, Yt) < si for some t < 2tiNNo. 

k<2NNQ 

We have shown that d{Xt, Yt) < e\ for some t < 2tiNNQ with probability greater 
thanp2 = Pi^^° > 0. By the Markov property, d(Xi, Yt) < e\ for some t < 2jtiNNQ with 
probability greater than 1 — (1 — P2)^ ■ To complete the proof, it suffices to let j — > 00. □ 



Two of the assumptions on the boundary of D made in Theorem 4.3, that it is convex 
with C-^-smooth boundary and it does not contain any line segments, are convenient from 
the technical point of view but most likely one can dispose of them with analysis more 
refined than that in our proof. 

Example 4.4. Suppose that D is the annulus {a; e : 1 < |a;| < 2}. The rotation of 
D around (0, 0) with an angle in (0, 27r) is an isometry with no fixed points. Hence, there 
exists a shy coupling of reflected Brownian motions in this annulus (see Example 3.3). 

There are many open problems concerning existence of shy couplings but we find the 
following two questions especially intriguing. Recall that B{x,r) denotes the open ball 
with center x and radius r. 
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Open problems 4.5. (i) Does there exist a shy couphng of reRected Brownian motions 
iiiB((0,0),3)\B((l,0),l)? 

(a) Does there exist a shy couphng of reflected Brownian motions in any simply 
connected planar domain? 

We end this paper with a vague remark concerning a potential relationship between 
shy couplings and an old and well known problem of "fixed points." Suppose that 5 is a 
topological space. If every continuous mapping I : S ^ S has a fixed point, i.e., a point 
X ^ S such that I{x) = x, then we say that S has the fixed point property. One of the 
most famous fixed point theorems is that of Brouwer — it asserts that a closed ball in M.'^ 
has the fixed point property. Spheres obviously do not have the fixed point property. Some 
of our results may suggest that a shy coupling exists if and only if the state space does 
not have the fixed point property. Example 3.8 applied to the graph illustrated in Fig. 3.6 
shows that this conjecture is false at this level generality. It is possible, though, that a 
weaker form of this assertion is true — we leave it as an open problem. 
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